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Persistence in systems with conserved order parameter 
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We consider the low-temperature coarsening dynamics of a one-dimensional Ising ferromagnet 
with conserved Kawasaki-like dynamics in the domain representation. Domains diffuse with size- 
dependent diffusion constant, D(l) oc P with 7 = — f. We generalize this model to arbitrary 7, 
and derive an expression for the domain density, N(t) ~ t - ^ with (f> = 1/(2 — 7), using a scaling 
argument. We also investigate numerically the persistence exponent 8 characterizing the power-law 
decay of the number, N p (t), of persistent (unflipped) spins at time t, and find N p (t) ~ t~ B where 9 
depends on 7. We show how the results for (f> and 6 are related to similar calculations in diffusion- 
limited cluster-cluster aggregation (DLCA) where clusters with size-dependent diffusion constant 
diffuse through an immobile 'empty' phase and aggregate irreversibly on impact. Simulations show 
that, while (j> is the same in both models, 6 is different except for 7 = 0. We also investigate models 
that interpolate between symmetric domain diffusion and DLCA. 

PACS numbers: 02.40.-r, 05.40.-a, 05.20.-y 
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I. INTRODUCTION 

Phase separation in the one-dimensional Ising ferro- 
magnet is a well- understood phenomenon Q, 0, At 
zero temperature, with Glauber (nonconserved) dynam- 
ics, coarsening occurs through the irreversible annihila- 
tion of domain walls, and the mean domain size grows 
as L(t) ~ t x l 2 . The coarsening process is equivalent to 
the diffusion-limited single-species annihilation reaction 
A+A — * 0, where A particles represent domain walls. Al- 
though no simple particle model can be constructed for 
Kawasaki (spin-exchange) dynamics due to the nearest- 
wall correlations introduced by the conservation princi- 
ple, qualitatively similar coarsening dynamics kinetics 
(but with a different domain-growth law, L(t) ~ t 1 / 3 ) 
have been observed for spin-exchange. 

Results have been obtained for the case of non- 
conserved Glauber dynamics for the growth law govern- 
ing the asymptotic increase of the average domain length 
L(t) and recently even an exact expression for the 'per- 
sistence' exponent 9 This exponent regulates the 
asymptotic decay, Q(t) ~ t~ e , of the persistence prob- 
ability, which is the probability that a spin has never 
flipped (or, equivalently, the fraction of spins which have 
not flipped) up to time t. It is a new dynamical exponent, 
not related to the domain-growth exponent Q. 

For conserved Kawasaki dynamics much less is known. 
It has been shown that the average domain length grows 
algebraically with time as L(t) ~ t 1 / 3 instead of ~ t 1 / 2 , 
a consequence of the l/l dependence that emerges at 
low temperature, from the underlying spin-exchange pro- 
cesses, for the effective diffusion constant of a domain of 
length I. Cornell and Bray p| have demonstrated that 
this size-dependence vanishes in the presence of an ex- 
ternal driving force. A calculation of the domain density, 
N(t), gave N(t) ~ < -1 / 2 as expected for size-independent 
dynamics. We will present a similar calculation for the 



unbiased Kawasaki model to obtain an expression for the 
density which we compare with data from numerical sim- 
ulations. 

To the best of our knowledge no analytical treatment of 
the persistence probability exists in the literature for the 
Kawasaki problem. The only numerical simulation of the 
system appears to be where the value of the exponent 
is measured as 9 — 0.73. We have repeated this study 
for a larger sample of initial configurations to obtain a 
similar value of 9, and have also simulated the same sys- 
tem with the general domain-diffusion law D(l) ~ P for 
a range of negative, integer 7. These models are not re- 
alisations of spin-exchange but are nonetheless perfectly 
legitimate in the domain description and are furthermore 
closely related to irreversible cluster-cluster aggregation 
(DLCA), another conserved system. We will show that 
our numerical values for the dynamic exponent <fi m this 
generalised Kawasaki model match the expression ob- 
tained by Helen et al. for DLCA, i.e. <\> cc = 1/(2-7), 
in the range that we examined. Interestingly, the bro- 
ken symmetry in cluster-cluster aggregation between the 
'filled' and 'empty' phases appears to have no effect on 
the rate of domain growth. It does, on the other hand, 
lead to differences between 9 CC — 2/(2 — 7) (where the 
superscript cc indicates cluster-cluster aggregation) and 
our own numerical data which yield, for 7 < —2, lower 
limits on #(7) in the symmetric spin system. 

The structure of the paper is as follows. In the next 
section we examine Kawasaki spin-exchange by reduc- 
ing the dynamics to domain diffusion in a suitable time- 
frame, and use a scaling hypothesis for the size distribu- 
tion of domains at time t, P(l,t), to calculate an expres- 
sion for the domain density in this regime. In section lll Bl 
we present data from simulations of this model which in- 
dicate that the persistence exponent 9 ~ 0.71. Sect ion IITTI 
introduces the other conserved system that will concern 
us, irreversible cluster-cluster aggregation, and contains 
a comparative treatment of both systems as well as sys- 
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terns with characteristics intermediate between these two 
models. We begin by introducing cluster-cluster aggre- 
gation and a few useful known results through which the 
correspondence between Kawasaki domain diffusion and 
aggregation cluster diffusion will be assessed, and con- 
tinue with an examination of the various exponents for 
weak (Sec lIII B|) and strong fSec lIII C|) size-dependence 
of the cluster (or domain) diffusion constant. 

II. THE KAWASAKI-ISING MODEL 

Unfortunately, a straightforward reduction to a par- 
ticle reaction-diffusion system, such as is the case for 
Glauber dynamics and the single-species annihilation 
process, is not possible for conserved dynamics. Exact 
mapping is retrieved in the infinite temperature limit, 
where the walls behave like branching-annihilating fixed 
particles [(| , but is rather complicated to follow at any fi- 
nite T. In the low-temperature limit of interest here, the 
simplest description of Kawasaki dynamics is in terms of 
domain diffusion 0, ITol | which we describe briefly below. 

The three spin processes that we wish to encapsulate 
in our higher-level description are: 

f. — | ► 1 — (spin diffusion) 

2. — | 1 ► h+ (spin capture) 

3. 1 — I > — I h (spin emission) 

where it is clear that capture is just the reverse of emis- 
sion and hence AE 3 = —AE 2 (and AEi = 0). The cor- 
responding rates differ by factors of order exp (—4 J/T) 
so, in the limit T <C J, the time-scales for diffusion and 
emission are widely separated, with capture occurring in- 
stantaneously. 

Implementing spin-exchange dynamics in this way on a 
ferromagnetic lattice from a random initial configuration 
leads initially to spin diffusion, followed by the capture 
of all available isolated spins until the latter disappear 
from the system. The rate of capture decays exponen- 
tially with time and contributes little towards coarsen- 
ing, whatever the size-distribution of domains at t = 0. 
Thereafter, the system becomes trapped in a metastable 
state whose lifetime is extremely long and diverges to 
infinity at zero temperature |l4j| . At any finite temper- 
ature, the metastability can only be broken once a rare 
splitting event occurs to provide a new isolated spin and 
even then the system refreezes quickly once the issued 
spin is absorbed. 

Spins produced in this way behave like random walk- 
ers, initially at i = 1, surrounded by two absorbing 
boundaries at x — and x — I, for a spin emitted into 
domain of length I. The spin is eventually captured by 
either of the boundaries with probabilities I — l/l and 1 jl 
respectively. Absorption at x — I results in the domain, 
across which the spin has diffused, being moved by one 
step whereas absorption at the other boundary returns 
the system to its original configuration. The combined 



rate for one step of the domain is ~ Z _1 exp (— 4 J/T) and 
it is this overall motion of the domain that facilitates 
coarsening close to T = Q • 

If t p is the physical time, in order to observe coarsening 
in the limit T — > 0, we must let t p — ► oo, while keeping 
i p exp(— AJ/T) fixed, and consider domain diffusion in 
rescaled time t = t p exp (—4 J/T). In this time variable, 
the rate of spin emission becomes of order unity, that 
of domain hops ~ l/l, where I is a typical domain size, 
and coarsening proceeds via domain annihilation induced 
by the diffusive dynamics of the domains. The complete 
set of microscopic spin processes has been reduced to a 
size-dependent domain diffusion model in time t. Cer- 
tain inconsistencies arise from this size-dependence for 
I = 1,2, but these 'end-effects' are negligible as far as 
the late-time dynamics is concerned (as is the early-time 
single-spin regime, which occurs instantaneously in terms 
of the time variable t) . Following the introduction of this 
coarse-grained model, we calculate the domain density 
N(t), and conclude the section with a discussion of the 
persistence probability Q(t). 

A. Domain density 

We focus first on the time evolution of the domain 
density by calculating N(t), the total number of domains 
per site at time t. The distribution of domain sizes, I, 
at time t is given by P(l,t), the number of domains of 
size I per lattice site. We expect that asymptotically the 
latter will be characterised by a single dynamical length 
scale, the average domain length L(t) ~ /r. The scaling 
hypothesis for P(l) is: 

P{l,t)~t-**P sc {x), x = ± (1) 

where the form of the prefactor is fixed by magnetisa- 
tion conservation, since J °° lP{l)dl must remain con- 
stant. The Z-dependence of the distribution at small 
I is a measure of the availability of domains near the 
I = boundary and dictates the rate of domain extinc- 
tion. More precisely, a domain disappears when it has 
length I = 1. Therefore dN/dt oc P(l,t) -> (dP/dl)\ i=0 
in the continuum limit. Since this rate must be fi- 
nite, it follows that P(l,t) oc I for I — > 0. A plot of 
Psc(x) = L 2 P(l,t) against x = l/L, where L is the 
mean domain size, confirms the linear small-argument 
behaviour (Fig^J whereas the same data plotted in the 
form ln(x _1 P sc ) versus x 2 (Fig|2J shows that the analytic 
form 

P sc {x) = Arexp(-A2 2 ) (2) 

provides a very good description of the data. 

It is straightforward to generalize the discussion to the 
case of non-zero magnetisation per spin Then the 
positive (+) and negative (— ) domains will have differ- 
ent distribution functions P±(l±,t). Numerical studies 
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FIG. 1: Plot of the scaling function P sc (a;) against the scaled 
length x from numerical simulations. Both the linear small- x 
and gaussian large- a; behaviour are reproduced. 




FIG. 2: The data of FigHlplotted instead as ln(P sc a; 1 ) versus 



confirm that the form J3J) provides a good fit to the data 
for + and — separately, but with different scale lengths 
L±(t). This suggests the scaling forms 



(3) 



where F(x) is defined such that F(0) = 0, and all nor- 
malisation factors have been combined in the constants 
a±. 

Returning to the rate of decay of N(t) we can write, in 



the general case, 



dN (. 1N dP + (l + ) 

, x . flp_(i_) 



(4) 



In the above expression we have used the fact that the 
rate for the elimination of + and — domains with I = 1 
is given by the mean rate at which — and + domains 
diffuse, i.e. (IZ 1 ) and respectively. The extra fac- 

tor of 2 comes from the fact that the elimination of a 
domain reduces the number of domains by 2 (since the 
domains neighbouring the eliminated domains merge). 
For an exact result, these averages should be taken only 
over domains adjacent to domains of size 1 = 1. In or- 
der to make further progress, we relax this condition and 
take the average over all domains of a given size. We are 
therefore neglecting correlations between domain sizes. 
The required averages (l^. 1 ) are given by 



(l±) 



J™l± 1 P±(l±)dl ± = 2a 
f™P±(l±)dl± 



N{t)L\ /o 



(5) 



where /„ = J °° x n F{x)dx and N(t) is the total number 
of domains per lattice site. The latter is obtained by 
adding the contributions N+(t) and N-(t) from + and 
— domains respectively: 



N(t) 
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/>Q<Q 

/ P_(l-)dl- = 
Jo 

2«+A 



la 
a_ 
1Z 



h 



(6) 



where the last equality follows from the condition 
N + (t) = N_(t) for all t, and holds for all values of 
the conserved quantity fi defined as the fraction of spins 
which have the value +1 (i.e. the magnetisation per site 
is 2fi — 1). This quantity satisfies 



and 



/>oo 

M = / l + p +( l +) dl + = a +h 
Jo 

poo 

1-H=l l_P_(l_)dl- = a-f 2 . 
Jo 



(7) 



(8) 



Using Eqs. an d © we can express (l^. 1 ), 

L± and a± in terms of N(t), \x and the integrals /„ (n — 
0, 1, 2). Putting these into Eq. (QJ and integrating yields 
the following expression for the domain density : 



N(t) 
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(9) 
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same for both domain types, and the neglect of correla- 
tions between domain sizes. 

To probe the quality of the fit more closely, we extract 
the amplitude C in the asymptotic behaviour N(t) — » 
Ct- 1 / 3 for large t. For fj, = 1/2, we find C = 0.440 ± 
0.003, whereas the analytic form gives C = 0.407, 
a 7.5% discrepancy. Ben-Nairn and Krapivsky [l(J have 
given an alternative treatment of the domain-size dis- 
tribution, making only the assumption that the domain 
sizes are uncorrelated. This leads to a slightly better 
value, 0.415, for C. However, no results were given for 
other values of [i in ref. O ur data show that the 

calculated C is smaller than the measured one by about 
7.5% for n = 0.3 and 0.1, and by about 5% for \i = 0.01. 
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FIG. 3: Time- dependence of the total domain density, N(t), 
for different volume fractions, fi, of the + phase. The con- 
tinuous curves show the analytical prediction given by Eqn. 
(reffinal). 

It remains to determine the constants /„. From the def- 
inition of the average domain lengths we find: 

J °° l±P±(l±)dl± _ h 

^ JT p±(i±)di ± fi 

=> h = fi ■ (io) 

For the assumed Gaussian scaling function F(x) = 
exp(— Ax 2 ), the required integrals give /o = (n/AA) 1 / 2 , 
fx = 1/(2A) and f 2 = (7T/16A 3 ) 1 / 2 . The result f 2 = fx 
fixes A = 7r/4, giving f = 1 and fx = fi = 2/tt. The 
final expression for N(t) becomes: 

N ® = UJfi-tf + nw) ' 

In Figure 01 we have plotted eq. (jl 1|) for different val- 
ues of /i together with numerical data from simulations 
described in the following subsection. The agreement is 
good even at early times, verifying the fast collapse to 
scaling observed in Figure ^ Note that in the deriva- 
tion of Eq. I|llf) we used the scaling form for all times 
t, not just late times. This may account for the small 
discrepancies between theory and data evident at early 
times in Figure [3] It is also evident that the data points 
lie slightly above the corresponding continuous curves at 
late times. Nevertheless, the generally good fit of Eq. 
(|llfl to the data, over a wide range of \i and with no 
adjustable parameters, provides good evidence that the 
assumptions made in its derivation are either correct or at 
least very good approximations. Principal among these 
assumptions are the Gaussian form taken for the scaling 
function F(x) (to derive explicit expressions for /o, fx 
and f-i), the assumption that this scaling function is the 



B. Persistence probability 

For the numerical simulations we adopted the domain 
description of the spin problem, since it is the most ap- 
propriate for T — > 0. We worked on ID lattices of size 
n = 10 5 with periodic boundary conditions. All sam- 
ples were prepared by randomly assigning spins to the n 
sites, although the system can be prepared in any other 
homogeneous initial state. The early-time form of the 
domain size-distribution is irrelevant to the asymptotic 
behaviour. 

Updating the persistence probability Q(t) means track- 
ing the fraction of spins that have not flipped or, equiv- 
alently, the fraction of sites that have not been visited 
by a domain wall up to t. On average, we expect each 
of the N domains to attempt to move by spin emission 
once per unit time in agreement with the rescaled rates, 
so a simple Monte-Carlo algorithm would be to select a 
domain at random, move it with probability 1 jl in either 
direction, where / is the its length, and increment time by 
At = 1/N(t) irrespective of whether or not the domain 
is moved. This is indeed the method used in Q to infer 
the value of the persistence exponent, 9 = 0.73. 

Unfortunately this algorithm suffers from severe 
slowing-down in the scaling limit, where the average do- 
main length L is much larger than unity and most hop- 
ping attempts are unsuccessful. We worked with a mod- 
ified version of the domain algorithm which does not 
exhibit such problems by ensuring that a move is ac- 
cepted at each computational step. A normalized vec- 
tor is constructed for the N domains, whose i th element, 
Hi(t) = ?i~ 1 /(Si=i T*)' ^ s tne probability that domain 
i experiences a hop at time t. For any given hop this 
probability is size-dependent only. Incrementing time by 
1/ "Y^—x l^ 1 1 the average time between successful moves, 
we end up with the same diffusion rate as in the sim- 
pler algorithm. We avoid, however, the problematic be- 
haviour relating to the inefficiency of the l/l trial. The 
disadvantage here is the calculation of the vector H at 
early times, which does result in the simpler algorithm 
being faster in this regime, but this is a small price to 
pay. We show in section ITlI CI that slow convergence to a 
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FIG. 4: Logarithm of the persistence probability against t' = 
Int. The local slope 9(t') is also given in the inset. 

scaling state in systems with diffusion rates D ~ P with 
7 < — 1 reinforces the need for the new algorithm. 

In Figure ^ we present a log-log plot of the persis- 
tence probability against time, averaged over 94 real- 
isations. The effective exponent 9(t) = —d\nQ/d\nt, 
plotted against lnt in the inset of Fig. refg-1, seems to 
saturate at value 9 ~ 0.716, the spread around this value 
beyond t ~ e 11 being primarily due to poor statistics 
resulting from the small number of persistent sites re- 
maining at this timescale. 

We note that, since 9 « 0.7 is larger than the domain- 
growth exponent (<p = 1/3), the typical distance between 
persistent sites, t e , grows faster than the domain length, 
t*. Clusters of non-persistent sites, of typical length t , 
grow by the motion of the domain walls at their bound- 
aries. These motions are uncorrelated since the separa- 
tion, t e , of the relevant walls is much larger than the do- 
main scale. The statistical independence of persistence- 
destroying processes suggests further that the distribu- 
tion function, rt(s), of the size s of non-persistence in- 
tervals should decay exponentially, n(s) oc exp(— s/(s)), 
where (s) ~ t e is the mean size of a non-persistent inter- 
val . This is to be contrasted with the Glauber model, 
in which the persistence length, i 3 / 8 , is smaller than the 
domain scale, i 1//2 , and n(s) has a nontrivial form [l5|. 

III. CONSERVED DYNAMICS IN 
AGGREGATION 

In this section we compare the coarsening dynamics 
of the Kawasaki model to that of aggregation models. 
How would the coarsening properties be modified in the 
Kawasaki ferromagnet for choices of the diffusion expo- 
nent 7 other than 7 = —1 in D ~ DqP7 Furthermore, is 
the amplitude ratio Rq = /Dq of the diffusion con- 



stants of the + and — domains relevant to the persistence 
probability and, if so, will a distinction between 'up' and 
'down' persistence be necessary? We discuss these points 
with the aid of results from irreversible cluster-cluster 
aggregation and extensive numerical simulations of mod- 
ified versions of the Kawasaki system. For the rest of 
the paper we shall use 'Kawasaki' to denote a more flex- 
ible model where we allow 7 / -1 and fij / 1. It 
must be pointed out, that both D ~ 1/Z and particularly 

Dq = Dq ' are robust features of Kawasaki dynamics, 
which follow directly from the elementary spin processes. 
This follows most strongly from the symmetric nature of 
the emission process where both an 'up' and a 'down' 
spin are emitted during a single event and consequently 
i?o = 1 is unavoidable. Although it is not possible to rec- 
oncile these rules with spin-exchange at the microscopic 
level, changes in the diffusion exponent and the ratio of 
amplitudes are easily accommodated in the domain de- 
scription that we have employed so far. 



A. Diffusion-limited cluster-cluster aggregation 
(DLCA) 

Irreversible cluster aggregation is closely related to 
Kawasaki dynamics. It, too, involves the separation of 
a two-phase mixture through domain diffusion and leads 
to conservation of its respective order parameter, in this 
case the total mass m. Clusters, defined as intervals of 
consecutive occupied sites, move randomly and coalesce 
irreversibly on impact, whilst the 'empty' phase through 
which the clusters diffuse is immobile (reversible aggre- 
gation scenarios, including various fragmentation mech- 
anisms or monomer 'chipping' Q, can actually result in 
effective diffusion of the empty phase, and hence partly 
restore the dynamic equivalence). Consequently, this 
model corresponds to the Rq = 'Kawasaki' system. 
The diffusion rate for a cluster of length I is Dp ~ I 1 
but we shall concentrate on 7 < and particularly on 
7 = 0,-1, which correspond to artificial but instructive 
realizations of domain diffusion on the Ising lattice. 

Due to the broken mobility symmetry between the 
'filled' (cluster) and 'empty' phases, a distinction should 
be made between the persistence properties of initially 
occupied sites and initially empty ones. Therefore, we 
use Pe(t) to denote the persistence probability of empty 
sites and Pp(i) that of occupied sites with similar nota- 
tion applying to the corresponding exponents. Only Pp 
is universal whereas Pp depends strongly on the volume 
fraction fi of the 'filled' phase, which of course remains 
constant throughout the coarsening process. 

For negative 7 the exponents governing the algebraic 
decay of the cluster density and empty persistence prob- 
ability have been calculated to be 4> cc = 1/(2 — 7) 
and 9e = 2</> cc for the DLCA model. The dependence 
of the dynamic exponent on 7 is readily obtainable from 
scaling alone |12| . an approach whose validity extends 
to 'Kawasaki' dynamics in view of sec III Al If the do- 
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FIG. 5: Log-log plot of the domain density with time for 
7 = —6 to —1 (top to bottom). The asymptotic slopes are 
-0.117,-0.136,-0.162,-0.197,-0.248 and -0.332 respectively, in 
reasonable agreement with the prediction <j> — 1/(2 — 7). 
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FIG. 7: Dynamic and persistence exponents for the size- 
independent 'Kawasaki' model in the range < Ro < 1. We 
have chosen De = D+. The upper and lower data sets give 
the persistence and growth exponents, 6 and (j>, respectively; 
6t refers to the persistence of the whole lattice and is obtained 
from the sum of the two constituent probabilities. The values 
of the exponents are consistent with <j> cc = 1/2 and 6 CC = 1. 




time required to reach the asymptotic value seems inde- 
pendent of the diffusion exponent 7 (FigJBJ, so that all 
0(7) saturate within the numerical window. The data 
suggest that domain growth is identical for the limiting 
cases Ro = (aggregation) and Ro = 1 ('Kawasaki'). We 
believe this to be the case in the whole range < Ro < 1 , 
given eq.Q and simulations indicating no systematic 
change in 0(7 = 0) for constant but unequal diffusion 
coefficients (Fig0). 
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FIG. 6: Plot of the difference (j>-cj> cc = (j)(Ro = l)-<f>(Ro = 0) 
in logarithmic time, for 3 values of 7. 



main scale is I, the time required to develop this scale 
is t ~ l 2 /D{l) ~ l 2 -~>, whence I ~ i 1 ^ 2 "^. For the 
empty-site persistence, i.e. the probability that an ini- 
tially empty site has not been covered by a cluster up to 
time t, Helen et al. used D(t) = Dot 1 * as a mean-field 
representation of the cluster kinetics to get 6e, which is 
known to be exact for 7 = 0. Numerical simulations indi- 
cate that the persistence of the 'filled' phase also decays 
algebraically with a non- universal exponent 0f(p) and 
that Of — 0e in the half- filled case, fi = 1/2, for all 7. 

Agreement with simulations of the zero-magnetisation 
'Kawasaki' model is very good for the form of the dy- 
namic exponent in the regime —1 < 7 < —6 (Fig[SJ). The 



B. The cases 7 = 0, -1 

Let us examine size- independent diffusion (7 = 0) first. 
The empty-site persistence probability for Rq — re- 
duces to the survival probability of a stationary point 
surrounded by two random walkers - the near edges of 
the two domains surrounding the point (FiglSJa)). This 
survival probability is given by the product of the sur- 
vival probability due to each walker and therefore decays 
asymptotically as t _1 / 2 < -1 / 2 = t~ l . Similarly, we can in- 
terpret filled-site persistence as the survival probability 
of a random walker with receding boundaries (Fig|8£b)). 
The persistent site carries the diffusion of the domain in 
which it lies, while the edges expand away from it due 
to domain collisions. For any other Ro the survival of 
the persistent site involves a mixture of both processes 
(FigEtc)). 

Although the exact kinetics of individual walls in 
FigJSJc) remain unknown, hints as to their form have 
emerged in many contexts. In a study of annihilating ran- 
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FIG. 8: Space-time diagram of persistence-destroying pro- 
cesses in (bottom to top) (a) empty, (b) filled and (c) sym- 
metric phases. 



dom walks on the infinite half-line it was shown from 
numerical measurements that the probability density of 
the outermost particle is asymmetric. In that case too 
the trajectory of the outermost walker is a combination of 
diffusion and directed jumps resulting from particle an- 
nihilation, similarly to Kawasaki walls. The probability 
density curve was found to decay like exp (—ax 2 ) in the 
direction of the free interface, but like exp (—(3x) on the 
side of the filled lattice, a behaviour which seems to sug- 
gest a regime where the diffusion exponential dominates 
and one where the annihilation exponential does. This 
appears reasonable on the basis that the contribution of 
these irregular hops is super-diffusive, characterised by 
the domain length distribution from which the jumps are 
drawn. In studies of the A + A — > reaction, this com- 
ponent has indeed been shown to exhibit an exponential 
form of the kind exp (— 8x)\ll\. Unfortunately, our cal- 
culation on a naive convolution of the two components 
has turned out an overestimate of the persistence expo- 
nent in the case where the persistent site is stationary 
and the probability can be expressed as the product of 
the two individual capture probabilities, suggesting that 
correlations between the diffusion and annihilation pro- 
cesses cannot be neglected. 

As far as numerical simulations of the size-independent 
system are concerned, we have found that within the ac- 
curacy of the statistics, no systematic shift in the values 
of the exponents takes place in the transition from ir- 
reversible aggregation to symmetric diffusion. It should 
be noted that the data in Fig[7| consistently fall slightly 
above the predicted value 9 — 1. Since the exponent is 
known to be 1 in the limiting cases Ro = 1 and Rq = 
Q it may be that the numerics have not yet reached 
the true asymptotic value. Nevertheless, even consider- 
ing this small discrepancy along all values of Rq, 7 = 



remains the only case for which the values of all four ex- 
ponents match in the limits Rq — and Rq = 1. For 
Kawasaki dynamics (7 = —1), which we discussed in 
Sec lII Bl a deviation of 9 from 9 CC is visible and begins 
to grow as we reduce the diffusion exponent further (see, 
Sec Iffl Cll . The inset in Figure 01 indicates that 9 sat- 
urates asymptotically above the 2/3 aggregation value, 
our calculations favouring 9 — 0.716(5). 

C. -2>7>-6 

In this parameter range, the aggregation and 
'Kawasaki' results for the persistence exponent clearly 
disagree. Our simulations did not reach, in most 
cases, the asymptotic regime (FigEJl. Nevertheless, they 
demonstrate that the exponents invariably exceed the 
aggregation prediction 8 CC — 2/(2 — 7). The very slow 
convergence we observe has, in fact, also arisen in other 
studies of aggregation- fragmentation models , a class 
to which the 'Kawasaki' model can be shown to relate. 
It was noted there, in numerical estimates of the cluster 
mass distribution function, that for 7 < —2 simulations 
failed to approach a scaling state in the time available. 
Our own efforts led us in many cases to final configu- 
rations where the average domain size was comparable 
to the lattice size. We think that for the algorithm em- 
ployed, good scaling-state statistics would require unfea- 
sibly long runs. We found no obvious way in which to 
improve further the performance of the routine which is 
still limited by the update of the probability vector H, al- 
though this penalty becomes less relevant at late times. 
Unfortunately, lattice renormalisation techniques (such 
as the one in [18J) which have proven successful in other 
circumstances can not be used in this case as they lead 
to local violations of the conservation principle. 

It is also clear that the collapse onto a single curve 
encountered in Figure at large t is absent for the per- 
sistence exponent. This much slower approach to asymp- 
totic behaviour for 9 than for (f> is also consistent with the 
statement made in Sec lII B| that non-persistent intervals 
and domain lengths have different scaling characteristics. 

Despite these complications, lower limits for #(7) are 
set by the data, and imply faster-than-0 cc decay for Q 
in the symmetric system. This is perhaps surprising in 
view of the results in the previous subsection, where such 
a deviation from the Rq = behaviour is absent. We 
do not know of any argument for the 7 = system to 
constitute a special case, but it is tempting to conjecture 
that the deviation of 9 from 9 CC observed for all 7 < —1 
is monotonic in the intermediate regime < Rq < 1. 

IV. DISCUSSION AND SUMMARY 

In this paper we considered an integrated analysis of 
two-phase conserved systems in one dimension based on 
two examples: irreversible cluster-cluster aggregation, 
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FIG. 9: Effective persistence exponent against lnt for Ro = 
1 'Kawasaki' systems with negative 7. The horizontal lines 
mark the values 6 CC predicted for Ro = 0. 



where only one of the phases is allowed to diffuse, and 
Kawasaki dynamics, where both phases diffuse at the 
same rate. In terms of the ratio of diffusion amplitudes 
Ro, the algebraic domain growth at late times has been 
established as identical in the these two limits, with the 
'Kawasaki' dynamic exponent (f> matching the prediction 
1/(2 — 7) whatever the form of size-dependence assumed 
for the diffusion rate. In the case of size-independent dif- 
fusion, we have also found that this equivalence holds for 
all other intermediate values of i?o, as was expected given 
the restrictions on scaling imposed by the conservation 
principle, i.e. that the order parameter should remain 
constant throughout coarsening. The robustness of this 
behaviour has been confirmed in a variety of aggregation 
systems 

By contrast, in the critical fi = 1/2 limit considered, 
changes in the relative mobility of the two phases lead to 
different persistence characteristics. The transition from 



irreversible aggregation to symmetric diffusion is accom- 
panied by severe slowing-down in the convergence rate of 
the effective persistence exponent to its asymptotic value 
as well as a distinct shift from the simple 9 CC = 2/(2 — 7) 
result for all negative 7. We have been able to show 
that the disagreement, which is at zero 7 absent in the 
whole range < Rq < 1, always leads to faster persis- 
tence decay for the symmetric system, with numerical 
data indicating a possible lower bound for the exponent 
at 9 = 2.30 (7 < -2). 

Apart from the slow convergence to asymptotic be- 
haviour which hinders any definitive results as far as the 
value of 9 in the symmetric limit is concerned, it is also 
possible that, for size-dependent cluster dynamics and 
finite values of Rq, Q(t) deviates from the simple power- 
law behaviour observed in irreversible aggregation. In 
this latter case, it seems that the filled persistence (the 
persistence of the mobile phase) is non-universal and de- 
cays as a power-law only for \i = 1/2 ||. Since single-site 
persistence for Ro > involves a mixture of both the 
empty- and filled-site behaviour, it is possible that the 
resulting probability is also non-universal. 

A further prospect which must be considered is the 
possibility of inferring 9 in the more strongly size- 
dependent systems from the scaling properties of the size- 
distribution of non-persistent sites. The general char- 
acteristics of this scaling are dictated by the fact that 
9 > <fi, which constrains the tail of the non-persistent 
interval size-distribution to a known form |15| . Hints of 
this behaviour have been suggested by our results for the 
probability Q(t), since the regularity encountered in the 
collapse of the growth exponent onto the aggregation re- 
sult for different 7 is lost in plots of 9. Instead, a very 
slow saturation of the local persistence exponent has been 
observed extending, in all cases but 7 = — 1, well beyond 
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